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1. SUMMARY 


The computer program described in this report was developed for use in 
a study of load sequence effects on crack growth, described in Reference l. 
Cycle by cycle stress histories were generated using techniques based on random 
noise theory, these techniques were implemented in a computer program described 
in Volume I, Reference 2. Using that program, stress histories were generated 
for three baseline spectra: Air-to-Air, Air-to-Ground, and Instrumentation 
and Navigation. Using the program described in this report, a combination of 
these was created to form the fourth baseline spectrum, the Composite. The 
stress histories for these four baseline spectra were modified to create 102 


spectra variations. 


2. PROGRAM OUTLINE 


This program combines and modifies stress history simulations, generated 
by a companion computer program, and creates four baseline spectra and 102 spectra 
variations. The baseline spectra are Air-to-Air (A-A), Air-to-Ground (A-G), 
Instrumentation and Navigation (I&N), and Composite. The spectra variation 
types are (a) Re-ordering of loads within a mission, (b) Sequence of missions, 
(c) Individual flight length, (d) Mission mix, (e) Individual flight length, 
(£) High and low load truncation, (g) Compression loads, (h) Exceedance curve 
variations, and (1) Coupling of peaks and valleys. Input to the program includes 
control numbers to create the spectra variations, and data sets of peak and valley 
stress history simulations obtained from the companion computer program, Volume I. 
Output includes a sequential list of peaks and valleys, and a table summarizing 
peak and valley coupling. Input and output are described in Sections 3 and 4, 
respectively. 

The program consists of a main program and eight subroutines, as outlined 


in Figure 1. The function of each sub-program is described in the following para- 


graphs. 


MAIN PROGRAM 


®@® Read input Contro! Data 
Print Input Control Data 


READER 
@ Reorder Stress in a @ Sum Peak and Valley @ Read Input List of 
Mission Hi-Lo or Lo-Hi Coupling Table Peaks and Valleys 
Perform Clipping 
Perform Truncation DISPLA 


increase Negative Stresses 
Add Lower Stresses 


@ Sum Columns and 
Rows in Peak and 
Valley Coupling Table 


@ Reorder Missions 


Hi-Lo or! Lo-Hi @ Print Peak and Valley 
Coupling Table 
Write Output List of Peaks. SEV 
HILOAD and Valleys 


@ Increase Peaks of 
@ Replace Stresses with Selected Missions by 
Higher Stresses a Factor 


FIGURE 1. SUBPROGRAMS USED IN MODIFYING STRESS TIME HISTORIES 


Main Program — The main program reads the input control data and prints 
the input data to permit a check. Subroutine READER is called to read the 
initial data sets of peak and valley stress history of each mission type (A-A, 

A~G, I&N) to be used in the spectrum. The main program is divided into three 
main process loops. Each loop develops the spectrum variation for a specific 
mission type either within the loop or by calling an appropriate subroutine. 

The three main process loops are very similar. They develop the flight- 
by-flight spectrum from the input peak and valley stress history by replacing 
selected valleys with stresses representing ground loads. As more peak and 
valley stress history is needed, subroutine reader is called and an additional 
data array of the required mission type is read. A series of IF tests, controlled 
by input data, determines the variations to be performed. The only variation 
performed within the main process loops is associated with coupling of peaks and 
valleys. The remainder of the variations are performed by subroutines SORTP, ORDER, 
SEV, HIL@AD and TAB. Messages at the end of the main program inform the user of 
input irregularities, such as a spectrum not ending with a stress representing 
ground load. 

Subroutine READER - This subroutine reads input sets of peak and valley 
stress history. The peak and valley sequence of stresses is maintained by deletion 
of the first stress value of the new data set, if necessary. The mission type is 
controlled by the input variable I@ (1 = A-A, 2 = A~G, 3 = I&N). 

Subroutine WRITEM - This subroutine is used to generate and write a sequential 
list of peaks and valleys. The array LIN contains a single line of data. IR is 
a count of number of values per line, with a maximum value of ten. 

Subroutine TAB - This subroutine has a dual purpose. First, it develops 
a table summarizing coupling of peaks and valleys which becomes a portion of program 
output. Secondly, it performs variations common to all mission types (A-A, A-G, 
I&N) such as clipping, increasing the magnitude of negative maneuver stresses, 
truncations and addition of low stresses. 

Subroutine DISPLA ~- This subroutine creates the column and line totals in the 
table summarizing coupling of peaks and valleys. 

Subroutine SEV - This subroutine increases the severity of the peaks in 
selected missions. 

Subroutine HILOAD - This subroutine replaces cycles with cycles containing 
the desired higher loads. 

Subroutine SORTP - This subroutine reorders the loads within each mission 


either in a Hi-Lo or Lo-Hi sequence. 


Subroutine ORDER - This subroutine reorders all missions into either a 


Hi-Lo or Lo-Hi sequence, based on the largest peak in each mission. 


3. INPUT DATA REQUIREMENTS 


Input to the program consists of control data to create the spectra 
variations, and three data sets containing peak and valley stress histories 
generated by companion computer program described in Volume I. These latter 
sets are identified as (1) A-A, (2) A-G, and (3) I&N. 

The first input control data are LINE] and GL@AD. These are used to des- 
cribe the mission mix (combination of A-A, A-G, and I&N) and, for each mission 
type, the stress representing the ground load. The completion of the mission 
mix description is signaled by input of the value END for LINE1. After comple- 
tion of the mission mix description, a series of forty-two indicators are read. 
Combinations of these indicators were used to create the 102 spectra variations 
studied in the work summarized in Reference l. 

The input data format is described in the computer program listing, via 
comment cards, Appendix A. Example input is presented in Appendix B; the data 
used to create the spectra studied in Reference 1 is listed in Appendix C. 


4. OUTPUT DESCRIPTION 

The output consists of three sets. The first set is an echo of the input 
control data to permit a check. 

The second set is a list of the cycle by cycle load spectrum; an option 
permits this spectrum to be written on a magnetic tape. The cycle by cycle 
load spectrum is listed both before and after reordering for spectrum varia~ 
tions with reordering of loads within a mission or with reordering of missions. 

The third data is a table summarizing the coupling of peaks and valleys. 
The table has the valley (% Design Limit Stress) as the ordinate and peak (4% 
Design Limit Stress) as the absicca. 


Example output is presented in Appendix B. 
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APPENDIX C 


INPUT DATA FOR STUDY SPECTRA VARIATIONS 
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APPENDIX D 


LIST OF COMPUTER PROGRAM SYMBOLS AND DEFINITIONS 
MAIN PROGRAM 


An integer counter for the number of items in the mission mix 
A do loop control variable used throughout the main program 
The main process loop control variable 

A do loop used when the multiple mission types are tallied 


A real array used to store one input array of Air-to-Air taw 
data 


A real array used to store one input array of Air-to-Ground 
raw data 


An integer mission mix counter 
A do loop control used when.printing the input data 


A real array used to store one input array of Instrument 
raw data 


A do loop control used to select the raw data from each raw 
data input array 


An integer line counter 


A do loop control used when an initial mission type does not 
start at the beginning of the raw data 


A do loop control used to read the four added peaks and valleys 
A do loop control used when printing the input data 

The number of load per mission 

An end of file check 

An escalation factor used to increase severity 

An integer peak and valley counter for Air-to-Air type missions 


An integer peak and valley counter for Air-to-Ground type 
missions 


A counter used to talley the seven possible added high loads 
An integer Instrument mission type peak and valley counter 


(See input parameters) 
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LIN 
LLL 
MIX 
SEV 
TAB 
TVP 
TVV 
ADDP 
ADDV 


FREQ 


IND1-IND9 


ISAA 


ISAG 


ISIN 


IV10-IV16 


IZIP 


NPAA 


NPAG 


NPIN 


ORDA 


PEAK 


VNM1 


AANM1 


PROGRAM SYMBOLS (MAIN) (Continued) 


The output line counter 

The mission type indicator 

The mission mix storage array 

The severity subroutine 

The occurrence table subroutine 

An occurrence table constraint for peaks 

An occurrence table constraint for valleys 

Used for adding additional high peaks 

Used for adding additional high valleys 

The interval used for adding additional high loads 

(See input parameters) 

A peak and valley counter switch for Air-to-Air type missions 
A peak and valley counter switch for Air-to-Ground type missions 
A peak and valley counter switch for Instrument type missions 
(See input parameters) 

An integer counter 

The number of points in an Air-to-Air input array 

The number of points in an Air-to-Ground input array 

The number of points in an Instrument input array 


A subroutine which orders or reorders missions by some preset 
convention 


The peaks used during a single mission 
V(N-1) the prior valley at the time of use 


A dummy variable used to insure the proper peak valley couple 
at the beginning of the first Air-to-Air input array 


A dummy variable used to insure the proper peak valley couple 
at the end of each Air-to-Ground input array 


AA(N-1) the last value of the prior Air-to-Air input array 
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AGDUM 


AGLAS 


AGNM1 


EOFAA 


EOFAG 


EOFIN 


GLOAD 


HWMNY 


INDAA 


INDAG 


INDIN 


INDUM 


IND1O-IND16 


INLAS 


INNM1 


INPAA 


INPAG 


INPIN 


IPICK 


LINE1 


LISTL 


PROGRAM SYMBOLS (MAIN) (Continued) 


A dummy variable used to insure the proper peak valley couple 
at the beginning of the first Air-to-Ground input array 


A dummy variable used to insure the proper peak valley couple 
at the end of each Air-to-Ground input array 


AG(N-1) the last value of the prior Air-to-Ground input array 
The end of file check for Air-to-Air missions 

The end of file check for Air-to-Ground missions 

The end of file check for Instrument missions 

The ground load 

A range number used in computing various mixes 

An integer switch used in selecting Air-to-Air input array 

An integer switch used in selecting Air-to-Ground input array 
An integer switch used in selecting Instrument input array 


A dummy variable used to insure the proper peak and valley 
couple at the beginning of the first Instrument input array 


(See input parameters) 


A dummy variable used to insure the proper peak and valley 
couple at the end of each Instrument input array 


IN(N-1) the last value of the prior Instrument input array 
For Air-to-Air type missions 

INPAA = (1) the next input array starts with a valley 
INPAA = (2) the next input array starts with a peak 

For Air-to-Ground type missions 

INPAG = (1) the next input array starts with a valley 
INPAG = (2) the next input array starts with a peak 

For Instrument type missions 

INPIN = (1) the next input array starts with a valley 
INPIN = (2) the next input array starts with a peak 


A binary switch used to select a peak or valley at any point 
in any mission 


An input carrier array 


This array is used to identify the input mission mix 
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LLINC 


MIXIT 


MTYPE 


ORDER 


SORTF 


SORTP 


TABLE 


WHERE 


DISPLA 


GLOADS 


HILOAD 


TERROR 


ISTART 


LLIBIT 


LODVAL 


MISION 


MISNUM 


READER 


VALLEY 


WRITEM 


PROGRAM SYMBOLS (MAIN) (Concluded) 


Counts the peaks and valleys in each mission 

Counts the types of missions in any mix 

Stores the mission types for use in the program 

A subroutine which reorders peaks and valleys within a mission 


An entry point into subroutine SORTP used when all missions 
have been entered and the sorting of missions ts about to occur. 


The mission ordering subroutine, entered at the end of each 
mission where the peaks and valleys of the mission are recorded 
on a direct access device for later sorting 


A subroutine which compiles the occurrence table and performs 
common mission variations 


The lower limit of each input array 


The subroutine with at the end of run prints the occurrence 
table 


The storage array for all ground loads 
The subroutine which adds selected high loads to missions 


A switch used to indicate that the final seed or mission does 
not end in a ground load 


An integer switch used when a desired mission does not start 
at the beginning of an input array 


Product of the number of cycles per mission and the number of 
missions 


The number of added high loads 

The mission counter 

The number of missions 

An input routine used to read and adjust each new input array 


A valley storage array used for reordering of peak valley pairs 
within a mission and for storing missions to be reordered 


A subroutine which is called each time a complete line of ten 
peaks and valleys are formed. 
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IO 


NP 


EOF 


IND 


LAS 


NM1 


VNM1 


PROGRAM SYMBOLS (READER) 


The A array contains one seed of peaks and valleys 


A control variable 

I = (1) the seed starts with a valley 
I = (2) the seed starts with a peak 
An integer constant 


The data set reference number 

IO = (1) is an Air-to-Air mission 

TO (2) is an Air-to-Ground mission 
IO (3) is an Instrument mission 


The number of data points in each seed 
The end of file indicator 


An indicator which shows the type of couple that was made 
between two seeds 


The last value of the new seed 
The next to last value of the new seed 


The last value of the prior seed 
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Ivl > 
LIN 
LLL 
TVP 
TVV 
VGL 
INDI > 
IV10 > 
ISWT 
GL@AD 
IND1O > 
NEXTV 


NPINC 


TABLE 


WRITEM 


PROGRAM SYMBOLS (TAB) 


A do loop control variable 

A do loop control variable 

The single peak that is entered via the call argument list 
The single valley that is entered via the call argument list 
An integer counter 

A counter which ten output values per line 

The number of points in a seed 

1V9 (See input parameters) 

Storage line for each print line 

An integer switch 

An occurrence table constraint for peak 

An occurrence table constraint for valleys 

No change in ground load check 

INDO (See input parameters) 

IV16 (See input parameters) 

An interger switch 

The ground load 

IND16 (See input parameters) 

The next logical valley 


The number of points in a seed times the frequency of points 
between ground loads 


The storage array for the occurrence table 


A subroutine which is called each time a complete line of ten 
peaks and valleys are formed 
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LLL 
DUM1 > 


MISION 


PROGRAM SYMBOLS (DISPLA) 


A do loop control variable 

A do loop control variable 

The storage array for the occurrence table 
Labels for the occurrence plot ordinate 
Labels for the occurrence plot absicca 

An integer subscript 

An integer subscript 

Dummy common value 

DUM3 dummy common value 


The mission counter 
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PROGRAM SYMBOLS (WRITEM) 


I A do loop control variable 

IO The data set reference number 

IR The number of data values per line 

LIN An array which contain ten or less output values 
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PL 

VL 
ADDP 
ADDV 
GLOAD 


LODVAL 


Pe a ee 


PROGRAM SYMBOLS (HILOAD) 


A counter used to limit the number of high loads added 


The 


The 


The 


The 


The 


The 


The 


peak for which the high load is substituted 
valley for which the high load is substituted 
line value of the peak 

line value of the valley 

added peak array 

added valley array 


ground load 


A same mission check 


A dummy common value not used in this subroutine 
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PROGRAM SYMBOLS (SEV) 


FAC The added severity factor 
PEAK The peak to be altered from the AA, AG or IN array 
PEAKL The peak to be altered from the print array LIN 


126 


IO 


IR 


IND 


IV1l 
IND1 
LIN 
IFLT 
GL@AD 
IPASS 
LLINC 
VALUE 


IFLTML 


PEAKMX 


PROGRAM SYMBOLS (ORDER) 


A do loop control used in search for the largest peak in each 
mission 


An array which holds the element number of the largest 
peak in each mission 


A do loop control used for printing the ordered mission 
A do loop control used in the sorting outer loop 

A do loop control used in the sorting inner loop 

The peaks for a single mission 

The valleys for a single mission 


The relative key used in reading the missions in an ordered 
manner 


The line item counter 


The part of the world is already in order indicated switch 
used in the sorting technique 


IV16 (See input parameters) 

IND16 (See input parameters) 

The print array used for printing each line of output 
The number of flights in the run 

The ground load 

Not used 

The number of peaks and valleys in a mission 

An intermediate value used in sorting 


The total number of flights minue one, used in the sorting 
technique 


An array which stores the maximum peak of each mission to be 
sorted at the conclusion of the run 
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IND 


Ivl + Iv9 
LIN 


LM1 


IND1 > IND16 
IV10 + IV16 
GLOAD 
LLINC 


TERROR 


PROGRAM SYMBOLS (SORTP) 


A do loop control variable used in the outer sorting loop 
A do loop control variable used in the inner sorting loop 
A sorting interchange temporary variable 


An integer array of size one hundred used in the sorting 
technique 


The peaks for a single mission 
The valleys for a single mission 
The line item counter 

Not used 


The part of the world is already in order indicator switch 
used in the sorting technique 


(See input parameters) 
The print array used for printing each line of input 


The number of values in a mission minus one used by the sorting 
technique 


(See input parameters) 

(See input parameters) 

The ground load 

The number of peaks and valleys in a mission 


The last value is a ground load valley indicator 
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